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We investigate the SU{2) Higgs model at a Higgs Boson mass of ~ 34 GeV for tem- 
peratures at the electroweak scale. We discuss in detail the critical temperature, 
the scalar field vacuum expectation value and the latent heat. We also consider 
for which temperatures the plasma can be regarded as radiation dominated. 



1 Introduction 

In the standard electroweak theory, the baryon violating processes are sup- 
pressed for temperatures below the electroweak scale (« 100 GeV). There- 
fore the currently observed baryon asymmetry was finally determined at this 
temperature scale. Since for ttiese temperatures infrared singularities render 
perturbation theory uncertainBcl, a non-perturbative treatment of the plasma 
is needed, which can be provided by numerical lattice simulations. For this 
purpose we performed large scale lattice simulations of the SU(2) Higgs model 
in 4 Euclidean dimensions. Our results also provide an estimation of the 
reliability LI of the reduction step, which is usedJxL relate simulation results 
obtained in effective models in three dimensions E3EiJ to physics. 

The results detailed in this talk are obtained with a Higgs boson mass 
of Mh — 34 GeV, for which a strong first order phase transition is observed 
at the electroweak scale. In the simulations we use the following action with 
inverse lattice spacings in the range 2Tc < < 5Tc 



S[U,v\ ^ /3^(l-iTrC/pO+E{5Tr(^^ 

pi X ^ 

+A [iTr {^t^x) - 1]' - Tr (^++^C/..^^.) . (1) 

Here Ux^fi denotes the gauge link variable, Upi the smallest Wilson loop and 
Lp^ is the scalar field in 2(g) 2 isospin matrix notation. 

2 Lines of constant physics and critical temperature 

To keep the renormalised couplings ga and Xr constant when stepping down in 
the lattice spacing a, we guess the bare /3 and A from the 1-loop renormalisation 
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Figure 1: Flow of the renormalised parameters with the lattice spacing. For Rhw '■= 
Mh/Mw we have QrRhw = (32Ah)^/^. The octagons refer to points with k shifted by 
5 • 10~^ resp. 4 ■ 10~^ in case of the finest lattice, when compared to the squares. 



group equation. The hopping parameter k is tuned to its critical value in 
a simulation at a temporal lattice extent Lt, as detailed in Ref. la. Which 
determines the critical temperature to be aTc = L^^ . 

In order to obtain the size of the lattice spacing in physical units, the 
Higgs and W-boson masses and the renormalised coupling 13 are measured at 
T = 0. In a detailed finite volume study at = 2Tc we find a spatial lattice 
extent Lg > 12, which relates to the correlation length of the Higgs boson 
as Ls > 2.9 ^H, to be sufficient to keep the renormalised parameters constant 
within an accuracy of 1% or even betterS. Using these physical volumes we 
investigate the infiuence of the finite a and the uncertainty of the critical 
hopping parameter on the renormalised parameters. The outcome is shown in 
figure H, where a small shift of k affects only the W-mass in lattice units. This 
leads to a significant contribution to the error of a from the uncertainty of 
K. The renormalised couplings and gfj however are unaffected within their 
small errors. A second important result is that the renormalised couplings 
are unchanged when the lattice spacing is decreased. This means that our 
simulations indeed follow the lines of constant physics. 

We are now in the position to give a precise estimate of the critical tem- 
perature in physical units. This is shown in figure The results from the 
previous paragraph are given by the white symbols. Their error bars contain 
the uncertainties of Kc and of both of the Mh values. Since the lattice artifacts 
of eq. ([l]) are 0{a^), we extrapolate with a quadratic ansatz to small a-values, 
which is justified from the good — 1- The result Tc/Mnla^o = 2.15(4) is 
shown by the filled symbol. Note the surprisingly small scaling violations of 
only 5%, when comparing the result at a^^ — 2Tc to the extrapolated value. 
With Mw ^ 80 GeV we get = 72.8 ± 1.3 GeV for a ^ 0. 

The perturbative estimate i shown by the dashed line, differs from the 
extrapolated lattice result by three standard deviations, so one cannot exclude 
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Figure 2: Extrapolation of Tc/Mh to the continuum limit wHch is given by the filled 
symbol. The dashed lines represents the perturbative predictions^. 



non- negligible higher order or non-perturbative contributions. 



3 Scalar field vacuum expectation value 



The scalar field vacuum expectation value v plays a prominent role in the 
semi-classical estimation of the sphaleron rate in the symmetry broken phase. 
A gauge invariant renormalised vacuum expectation value is defined by := 
2Mw/gR- However since the mass determination is demanding in terms of 
CPU time, we approximate this by v{T) [2k {pI{T) - pl,^y,,^{Tc))Y/^ with 
:= iTr((/3+(^j;). It has been ijhecked that vr and v{T) are in reasonable 
agreement for T = and T k, TcU. The result for v /T is given in figure |3[ We 
also given an estimate of the supercooling obtained in the thin wall approxi- 
mation. From this one can expect at most an increase of 15% in the exponent 
of the sphaleron rate due to supercooling. With the Clausius-Clapeyron Equa- 
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tioni the latent heat can be estimated from the discontinuity at Tj,- 

Ae/T^ = 0.281(19)U-i^2T. , Ae/T^ = 0.31(12) U-i=4T. • (2) 
4 Thermodynamics of the plasma 

In this section we discuss the thermodynamic quantity S := — P which 
measures the deviation of the plasma from pure radiation. With e the energy 
density and with P the pressure is denoted, while S can be determined from 



S 



(3) 



with T := — log(Mvi/) and 



:= ^Tr ((/?^^^[/^_^(/5j.). We determine the deriva- 



tives of P and A using the 1-loop renormalisation group equations and those 
of K from fits to the simulation results for Kc u. Since eq. (^) contains diver- 
gent vacuum contributions, S{T — 0) must be subtracted to obtain the physical 
result, which is show in figure ^. The figure shows different errors: statistical 
errors by vertical bars, uncertainty of Kc by horizontal bars and the uncertainty 
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Figure 4: Result for (^e — P)T~^. Triangles represent = 2Tc, squares 4Tc. 

arising from ^ by the shaded region. In the symmetric phase 6/T goes down 
by one order of magnitude when T is increased from Tc to 2Tc where it is al- 
most compatible with 0. This strong jdecay is confirmed by the good agreement 
of the slopes for both a- values B, so that for T > 2Tc the plasma can be 
considered to be radiation dominated. 

Again from the jump at Tc the latent heat can be determined: 

Ae/T,4 = 0.240(30 + 4) U-i=2T., Ae/T^ ^ 0.28(3 + 9) U-i^4T. • (4) 

The first number in the brackets represents the statistical error and the second 
one the uncertainty of Kc- The result is in agreement to eq. (§) and no scaling 
violation is to be observed between the results at = 2Tc and ATc- 
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